Protection from successive Omicron variants with SARS-CoV-2 vaccine and monoclonal antibodies in kidney transplant recipients.

Introduction: Kidney transplant recipients (KTRs) are at high risk of severe COVID-19, even when they are fully vaccinated. Additional booster vaccinations or passive immunization with prophylactic monoclonal antibodies are recommended to increase their protection against severe COVID-19. Methods: Here, we describe the neutralization of SARS-CoV-2 Delta, Omicron BA.1, BA.2, BA.4, and BA.5 variants, firstly by 39 serum samples from vaccinated KTRs exhibiting anti-spike antibody concentrations ≥264 binding antibody units (BAU)/mL and, secondly, by tixagevimab/cilgavimab. Results: No neutralization was observed for 18% of the KTRs, while serum from only 46% of patients could neutralize the five variants. Cross-neutralization of the Delta and Omicron variants occurred for 65-87% of sera samples. The anti-spike antibody concentration correlated with neutralization activity for all the variants. The neutralization titers against the Delta variant were higher in vaccinated KTRs who had previously presented with COVID-19, compared to those KTRs who had only been vaccinated. Breakthrough infections occurred in 39% of the KTRs after the study. Tixagevimab/cilgavimab poorly neutralizes Omicron variants, particularly BA.5, and does not neutralize BQ.1, which is currently the most prevalent strain. Discussion: As a result, sera from seropositive vaccinated KTRs had poor neutralization of the successive Omicron variants. Several Omicron variants are able to escape tixagevimab/cilgavimab.

Correlation of Redondovirus and Entamoeba gingivalis Detections in the Human Oral Cavity Suggests That This Amoeba Is Possibly the Redondovirus Host.

The virome of the human oral cavity and the relationships between viruses and diseases such as periodontitis are scarcely deciphered. Redondoviruses were reported in the human oral cavity in 2019, including in periodontitis patients. Here, we aimed at detecting redondoviruses and at searching for a potential viral host in human saliva. Non-stimulated saliva was collected between December 2020 and June 2021. These samples were tested using real-time PCR regarding the presence of redondovirus and Entamoeba gingivalis DNA. Similarity searches were performed using BLAST against eukaryotic and prokaryotic sequences from GenBank. The redondovirus DNA was detected in 46% of the 28 human saliva samples. In addition, short fragments of redondovirus genomes were detected in silico within Entamoeba sequences. Finally, Entamoeba gingivalis DNA was detected in 46% of the 28 saliva samples, with a strong correlation between redondovirus DNA and E. gingivalis DNA detections, in 93% of the cases. Regarded together, these findings and previous ones strongly support the presence of redondoviruses in the human oral cavity and their association to E. gingivalis as their likely host.

High Throughput SARS-CoV-2 Genome Sequencing from 384 Respiratory Samples Using the Illumina COVIDSeq Protocol.

The emergence of the Coronavirus Disease 2019 (COVID-19) pandemic has fostered the use of high-throughput techniques to sequence the entire severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome and track its evolution. The present study proposes a rapid and relatively less expensive sequencing protocol for 384 samples by adapting the use of an Illumina NovaSeq library to an Illumina MiSeq flow cell instrument. The SARS-CoV-2 genome sequences obtained with Illumina NovaSeq and those obtained using MiSeq instruments were compared with the objective to validate the new, modified protocol. A total of 356 (94.6%) samples yielded interpretable sequences using the modified Illumina COVIDSeq protocol, with an average coverage of 91.6%. By comparison, 357 (94.9%) samples yielded interpretable sequences with the standard COVIDSeq protocol, with an average coverage of 95.6%. Our modified COVIDSeq protocol could save 14,155 euros per run and yield results from 384 samples in 53.5 h, compared to four times 55.5 h with the standard Illumina MiSeq protocol. The modified COVIDSeq protocol thus provides high quality results comparable to those obtained with the standard COVIDSeq protocol, four times faster, while saving money.

The emergence, spread and vanishing of a French SARS-CoV-2 variant exemplifies the fate of RNA virus epidemics and obeys the Mistigri rule.

The nature and dynamics of mutations associated with the emergence, spread, and vanishing of SARS-CoV-2 variants causing successive waves are complex. We determined the kinetics of the most common French variant ("Marseille-4") for 10 months since its onset in July 2020. Here, we analyzed and classified into subvariants and lineages 7453 genomes obtained by next-generation sequencing. We identified two subvariants, Marseille-4A, which contains 22 different lineages of at least 50 genomes, and Marseille-4B. Their average lifetime was 4.1 ± 1.4 months, during which 4.1 ± 2.6 mutations accumulated. Growth rate was 0.079 ± 0.045, varying from 0.010 to 0.173. Most of the lineages exhibited a bell-shaped distribution. Several beneficial mutations at unpredicted sites initiated a new outbreak, while the accumulation of other mutations resulted in more viral heterogenicity, increased diversity and vanishing of the lineages. Marseille-4B emerged when the other Marseille-4 lineages vanished. Its ORF8 gene was knocked out by a stop codon, as reported in SARS-CoV-2 of mink and in the Alpha variant. This subvariant was associated with increased hospitalization and death rates, suggesting that ORF8 is a nonvirulence gene. We speculate that the observed heterogenicity of a lineage may predict the end of the outbreak.

Simple Wastewater Preparation Protocol Applied to Monitor the Emergence of the Omicron 21L/BA.2 Variant by Genome Sequencing.

Detecting and monitoring viruses in wastewater samples have been reported as useful ways of tracking SARS-CoV-2 epidemic trends. However, there is currently no unanimously recognised method of processing samples to identify and quantify SARS-CoV-2 variants in wastewater. We aimed to implement a method that was as simple as possible in order to be used universally. In a study performed between January 2022 and June 2022 in the city of Marseille, France, we first evaluated the impact of the sample preservation strategy. We then compared ultracentrifugation to ultrafiltration and several steps of filtration to determine the optimal approach for virus concentration. As a proof-of-concept, the definitive protocol was applied to next-generation sequencing of SARS-CoV-2 in wastewater to monitor the emergence of the Omicron variant in the city. For sewage water to be processed in the week following the sampling, storage at +4 °C is sufficient, with less than 1 Ct loss. Filtration with a 5 µm syringe filter, then with a 0.8 µm filtration unit, followed by ultrafiltration was the optimal protocol, leading to an average increase of 3.24 Ct when the starting Ct was on average 38 in the wastewater. This made it possible to observe the emergence of the Omicron 21L/BA.2 variant after Omicron 21K/BA.1 by genome sequencing over a period ranging from 20 February to 10 April 2022 in agreement with observations based on patient data. To conclude, by using a simple method requiring only basic filters and a centrifuge as equipment, it is possible to accurately track the relative incidence rates and the emergence of SARS-CoV-2 variants based on sewage samples.

Quasispecies Analysis of SARS-CoV-2 of 15 Different Lineages during the First Year of the Pandemic Prompts Scratching under the Surface of Consensus Genome Sequences.

The tremendous majority of SARS-CoV-2 genomic data so far neglected intra-host genetic diversity. Here, we studied SARS-CoV-2 quasispecies based on data generated by next-generation sequencing (NGS) of complete genomes. SARS-CoV-2 raw NGS data had been generated for nasopharyngeal samples collected between March 2020 and February 2021 by the Illumina technology on a MiSeq instrument, without prior PCR amplification. To analyze viral quasispecies, we designed and implemented an in-house Excel file ("QuasiS") that can characterize intra-sample nucleotide diversity along the genomes using data of the mapping of NGS reads. We compared intra-sample genetic diversity and global genetic diversity available from Nextstrain. Hierarchical clustering of all samples based on the intra-sample genetic diversity was performed and visualized with the Morpheus web application. NGS mapping data from 110 SARS-CoV-2-positive respiratory samples characterized by a mean depth of 169 NGS reads/nucleotide position and for which consensus genomes that had been obtained were classified into 15 viral lineages were analyzed. Mean intra-sample nucleotide diversity was 0.21 ± 0.65%, and 5357 positions (17.9%) exhibited significant (>4%) diversity, in ≥2 genomes for 1730 (5.8%) of them. ORF10, spike, and N genes had the highest number of positions exhibiting diversity (0.56%, 0.34%, and 0.24%, respectively). Nine hot spots of intra-sample diversity were identified in the SARS-CoV-2 NSP6, NSP12, ORF8, and N genes. Hierarchical clustering delineated a set of six genomes of different lineages characterized by 920 positions exhibiting intra-sample diversity. In addition, 118 nucleotide positions (0.4%) exhibited diversity at both intra- and inter-patient levels. Overall, the present study illustrates that the SARS-CoV-2 consensus genome sequences are only an incomplete and imperfect representation of the entire viral population infecting a patient, and that quasispecies analysis may allow deciphering more accurately the viral evolutionary pathways.

SARS-CoV-2 Spike Protein Induces Hemagglutination: Implications for COVID-19 Morbidities and Therapeutics and for Vaccine Adverse Effects.

Experimental findings for SARS-CoV-2 related to the glycan biochemistry of coronaviruses indicate that attachments from spike protein to glycoconjugates on the surfaces of red blood cells (RBCs), other blood cells and endothelial cells are key to the infectivity and morbidity of COVID-19. To provide further insight into these glycan attachments and their potential clinical relevance, the classic hemagglutination (HA) assay was applied using spike protein from the Wuhan, Alpha, Delta and Omicron B.1.1.529 lineages of SARS-CoV-2 mixed with human RBCs. The electrostatic potential of the central region of spike protein from these four lineages was studied through molecular modeling simulations. Inhibition of spike protein-induced HA was tested using the macrocyclic lactone ivermectin (IVM), which is indicated to bind strongly to SARS-CoV-2 spike protein glycan sites. The results of these experiments were, first, that spike protein from these four lineages of SARS-CoV-2 induced HA. Omicron induced HA at a significantly lower threshold concentration of spike protein than the three prior lineages and was much more electropositive on its central spike protein region. IVM blocked HA when added to RBCs prior to spike protein and reversed HA when added afterward. These results validate and extend prior findings on the role of glycan bindings of viral spike protein in COVID-19. They furthermore suggest therapeutic options using competitive glycan-binding agents such as IVM and may help elucidate rare serious adverse effects (AEs) associated with COVID-19 mRNA vaccines, which use spike protein as the generated antigen.

Matrix-Encoding Gene Diversity of 624 Influenza A/H3N2 Genomes Does Not Show Association with Impaired Viral Detection by Commercialized qPCR Assays.

As for the case of SARS-CoV-2, genome sequencing of influenza viruses is of potential interest to raise and address virological issues. Recently, false-negativity of real-time reverse transcription-PCR (qPCR) assays that detect influenza A/H3N2 virus RNA were reported and associated with two mutations (A37T and C161T) in the Matrix-encoding (M1) gene located on viral segment 7. This triggered a national alert in France. The present study sought to assess the association between the presence of these mutations and potential false negative results of influenza A/H3N2 virus RNA detection by commercialized qPCR assays at the clinical virology laboratory of our university hospitals in southern France. This study focused on the genetic diversity in the M1 gene and segment 7 of 624 influenza A/H3N2 virus genomes obtained from respiratory samples having tested qPCR-positive with M1 gene-targeting assays in our clinical virology laboratory. A total of 585 among the 624 influenza A/H3N2 virus genomes (93.7%) were of clade 3C.2a1b.2a.2, and 39 (6.3%) were of clade 3C.2a1b.1a. M1 gene substitutions A37T and C161T were both present in 582 (93.3%) genomes, only of clade 3C.2a1b.2a.2. Substitution A37T was present in 621 (99.5%) genomes. Substitution C161T was present in 585 genomes (93.8%), all of clade 3C.2a1b.2a.2. Moreover, 21 other nucleotide positions were mutated in ≥90% of the genomes. The present study shows that A37T/C and C161T mutations, and other mutations in the M1 gene and segment 7, were widely present in influenza A/H3N2 virus genomes recovered from respiratory samples diagnosed qPCR-positive with commercialized assays.

First cases of infection with the 21L/BA.2 Omicron variant in Marseille, France.

The SARS-CoV-2 21K/BA.1, 21L/BA.2, and BA.3 Omicron variants have recently emerged worldwide. To date, the 21L/BA.2 Omicron variant has remained very minority globally but became predominant in Denmark instead of the 21K/BA.1 variant. Here, we describe the first cases diagnosed with this variant in south-eastern France. We identified 13 cases using variant-specific qPCR and next-generation sequencing between 28/11/2021 and 31/01/2022, the first two cases being diagnosed in travelers returning from Tanzania. Overall, viral genomes displayed a mean (±standard deviation) number of 65.9 ± 2.5 (range, 61-69) nucleotide substitutions and 31.0 ± 8.3 (27-50) nucleotide deletions, resulting in 49.6 ± 2.2 (45-52) amino acid substitutions (including 28 in the spike protein) and 12.4 ± 1.1 (12-15) amino acid deletions. Phylogeny showed the distribution in three different clusters of these genomes, which were most closely related to genomes from England and South Africa, from Singapore and Nepal, or from France and Denmark. Structural predictions highlighted a significant enlargement and flattening of the surface of the 21L/BA.2 N-terminal domain of the spike protein compared to that of the 21K/BA.1 Omicron variant, which may facilitate initial viral interactions with lipid rafts. Close surveillance is needed at global, country, and center scales to monitor the incidence and clinical outcome of the 21L/BA.2 Omicron variant.

Choosing a cellular model to study SARS-CoV-2.

As new pathogens emerge, new challenges must be faced. This is no different in infectious disease research, where identifying the best tools available in laboratories to conduct an investigation can, at least initially, be particularly complicated. However, in the context of an emerging virus, such as SARS-CoV-2, which was recently detected in China and has become a global threat to healthcare systems, developing models of infection and pathogenesis is urgently required. Cell-based approaches are crucial to understanding coronavirus infection biology, growth kinetics, and tropism. Usually, laboratory cell lines are the first line in experimental models to study viral pathogenicity and perform assays aimed at screening antiviral compounds which are efficient at blocking the replication of emerging viruses, saving time and resources, reducing the use of experimental animals. However, determining the ideal cell type can be challenging, especially when several researchers have to adapt their studies to specific requirements. This review strives to guide scientists who are venturing into studying SARS-CoV-2 and help them choose the right cellular models. It revisits basic concepts of virology and presents the currently available in vitro models, their advantages and disadvantages, and the known consequences of each choice.

SARS-CoV-2 quasi-species analysis from patients with persistent nasopharyngeal shedding.

At the time of a new and unprecedented viral pandemic, many questions are being asked about the genomic evolution of SARS-CoV-2 and the emergence of different variants, leading to therapeutic and immune evasion and survival of this genetically highly labile RNA virus. The nasopharyngeal persistence of infectious virus beyond 17 days proves its constant interaction with the human immune system and increases the intra-individual mutational possibilities. We performed a prospective high-throughput sequencing study (ARTIC Nanopore) of SARS-CoV-2 from so-called "persistent" patients, comparing them with a non-persistent population, and analyzing the quasi-species present in a single sample at time t. Global intra-individual variability in persistent patients was found to be higher than in controls (mean 5.3%, Standard deviation 0.9 versus 4.6% SD 0.3, respectively, p < 0.001). In the detailed analysis, we found a greater difference between persistent and non-persistent patients with non-severe COVID 19, and between the two groups infected with clade 20A. Furthermore, we found minority N501Y and P681H mutation clouds in all patients, with no significant differences found both groups. The question of the SARS-CoV-2 viral variants' genesis remains to be further investigated, with the need to prevent new viral propagations and their consequences, and quasi-species analysis could be an important key to watch out.

Limited permissibility of ENL-R and Mv-1-Lu mink cell lines to SARS-CoV-2.

The SARS-CoV-2 pandemic started in the end of 2019 in Wuhan, China, which highlighted the scenario of frequent cross-species transmission events. From the outbreak possibly initiated by viral spill-over into humans from an animal reservoir, now we face the human host moving globally while interacting with domesticated and peridomestic animals. The emergence of a new virus into the ecosystem leads to selecting forces and species-specific adaptations. The adaptation of SARS-CoV-2 to other animals represents a risk to controlling the dissemination of this coronavirus and the emergence of new variants. Since 2020, several mink farms in Europe and the United States have had SARS-CoV-2 outbreaks with human-mink and mink-human transmission, where the mink-selected variants possibly hold evolutionary concerning advantages. Here we investigated the permissibility of mink lung-derived cells using two cell lines, Mv-1-Lu and ENL-R, against several lineages of SARS-CoV-2, including some classified as variants of concern. The viral release rate and the infectious titers indicate that these cells support infections by different SARS-CoV-2 lineages. The viral production occurs in the first few days after infection with the low viral release by these mink cells, which is often absent for the omicron variant for lung cells. The electron microscopy reveals that during the viral replication cycle, the endomembrane system of the mink-host cell undergoes typical changes while the viral particles are produced, especially in the first days of infection. Therefore, even if limited, mink lung cells may represent a selecting source for SARS-CoV-2 variants, impacting their transmissibility and pathogenicity and making it difficult to control this new coronavirus.

First Detection of the SARS-CoV-2 Omicron BA.5/22B in Monaco.

The Omicron BA.5/22B variant has been designated as a "variant of concern" by the World Health Organization. We describe, here, the first evidence in Monaco of infection with an Omicron BA.5/22B variant, probably imported from the Republic of Seychelles, harboring a rare combination of non-BA.5/22B signature amino acid changes. SARS-CoV-2 neutralizing antibodies were measured with a surrogate virus neutralization test. SARS-CoV-2 genotype screening was performed on nasopharyngeal samples with a multiplex qPCR assay. The SARS-CoV-2 genome was obtained by next-generation sequencing with the Illumina COVID-seq protocol, then assembly using bioinformatics pipelines and software was performed. The BA.5/22B spike protein structure was obtained by molecular modeling. Two spouses were SARS-CoV-2-diagnosed the day they returned from a one-week trip in the Republic of Seychelles. SARS-CoV-2 qPCR screening for variant-specific mutations identified an Omicron variant BA.1/21K, BA.4/22A, or BA.5/22B. A SARS-Co-2 BA.5/22B variant genome was recovered from one of the spouses. Aside from BA.5/22B-defining amino acid substitutions, four other amino acid changes were encoded including Q556K in ORF1a, K2557R in ORF1b, and A67V and A829T in spike; only 13 genomes in sequence databases harbored these four mutations concurrently. Structural analysis of this BA.5/22B variant predicted that A829T in spike may result in a compaction that may affect conformational plasticity. Overall, our findings warrant performing genome-based genotypic surveillance to survey accurately the emergence and circulation of SARS-CoV-2 variants worldwide and point out that their first occurrence in a country is often through international travel despite implemented countermeasures.

Reinfections with Different SARS-CoV-2 Omicron Subvariants, France.

We describe 188 patients in France who were successively infected with different SARS-CoV-2 Omicron subvariants, including BA.1, BA.2, and BA.5. Time between 2 infections was <90 days for 50 (26.6%) patients and <60 days for 28 (14.9%) patients. This finding suggests that definitions for SARS-CoV-2 reinfection require revision.

Simultaneous SARS-CoV-2 Genome Sequencing of 384 Samples on an Illumina MiSeq Instrument through Protocol Optimization.

In the present study, we propose a high-throughput sequencing protocol using aNextera XT Library DNA kit on an Illumina MiSeq instrument. We made major modifications to this library preparation in order to multiplex 384 samples in a single Illumina flow cell. To validate our protocol, we compared the sequences obtained with the modified Illumina protocol to those obtained with the GridION Nanopore protocol. For the modified Illumina protocol, our results showed that 94.9% (357/376) of the sequences were interpretable, with a viral genome coverage between 50.5% and 99.9% and an average depth of 421×. For the GridION Nanopore protocol, 94.6% (356/376) of the sequences were interpretable, with a viral genome coverage between 7.0% and 98.6% and an average depth of 2123×. The modified Illumina protocol allows for gaining EUR 4744 and returning results of 384 samples in 53.5 h versus four times 55.5 h with the standard Illumina protocol. Our modified MiSeq protocol yields similar genome sequence data as the GridION Nanopore protocol and has the advantage of being able to handle four times more samples simultaneously and hence is much less expensive.

A 21L/BA.2-21K/BA.1 "MixOmicron" SARS-CoV-2 hybrid undetected by qPCR that screen for variant in routine diagnosis.

Among the multiple SARS-CoV-2 variants identified since summer 2020, several have co-circulated, creating opportunities for coinfections and potentially genetic recombinations that are common in coronaviruses. Viral recombinants are indeed beginning to be reported more frequently. Here, we describe a new SARS-CoV-2 recombinant genome that is mostly that of a Omicron 21L/BA.2 variant but with a 3' tip originating from a Omicron 21K/BA.1 variant. Two such genomes were obtained in our institute from adults sampled in February 2022 in university hospitals of Marseille, southern France, by next-generation sequencing carried out with the Illumina or Nanopore technologies. The recombination site was located between nucleotides 26,858-27,382. In the two genomic assemblies, mean sequencing depth at mutation-harboring positions was 271 and 1362 reads and mean prevalence of the majoritary nucleotide was 99.3 ± 2.2% and 98.8 ± 1.6%, respectively. Phylogeny generated trees with slightly different topologies according to whether genomes analyzed were depleted or not of the 3' tip. This 3' terminal end brought in the Omicron 21L/BA.2 genome a short transposable element of 41 nucleotides named S2m that is present in most SARS-CoV-2 except a few variants among which the Omicron 21L/BA.2 variant and may be involved in virulence. Importantly, this recombinant is not detected by currently used qPCR that screen for variants in routine diagnosis. The present observation emphasizes the need to survey closely the genetic pathways of SARS-CoV-2 variability by whole genome sequencing, and it could contribute to gain a better understanding of factors that lead to observed differences between epidemic potentials of the different variants.

The Role of Lebanon in the COVID-19 Butterfly Effect: The B.1.398 Example.

In the present study, we provide a retrospective genomic surveillance of the SARS-CoV-2 pandemic in Lebanon; we newly sequence the viral genomes of 200 nasopharyngeal samples collected between July 2020 and February 2021 from patients in different regions of Lebanon and from travelers crossing the Lebanese-Syrian border, and we also analyze the Lebanese genomic dataset available at GISAID. Our results show that SARS-CoV-2 infections in Lebanon during this period were shaped by the turnovers of four dominant SARS-CoV-2 lineages, with B.1.398 being the first to thoroughly dominate. Lebanon acted as a dispersal center of B.1.398 to other countries, with intercontinental transmissions being more common than within-continent. Within the country, the district of Tripoli, which was the source of 43% of the total B.1.398 sequences in our study, was identified as being an important source of dispersal in the country. In conclusion, our findings exemplify the butterfly effect, by which a lineage that emerges in a small area can be spread around the world, and highlight the potential role of developing countries in the emergence of new variants.

SARS-CoV-2 Testing of Aircraft Wastewater Shows That Mandatory Tests and Vaccination Pass before Boarding Did Not Prevent Massive Importation of Omicron Variant into Europe.

BACKGROUND: Most new SARS-CoV-2 epidemics in France occurred following the importation from abroad of emerging viral variants. Currently, the risk of new variants being imported is controlled based on a negative screening test (PCR or antigenic) and proof of up-to-date vaccine status, such as the International Air Transport Association travel pass. METHODS: The wastewater from two planes arriving in Marseille (France) from Addis Ababa (Ethiopia) in December 2021 was tested by RT-PCR to detect SARS-CoV2 and screen for variants. These tests were carried out between landing and customs clearance and were then sequenced by MiSeq Illumina. Antigenic tests and sequencing by NovaSeq were carried out on respiratory samples collected from the 56 passengers on the second flight. RESULTS: SARS-CoV-2 RNA suspected of being from the Omicron BA.1 variant was detected in the aircraft's wastewater. SARS-CoV2 RNA was detected in 11 [20%) passengers and the Omicron BA.1 variant was identified. CONCLUSION: Our work shows the efficiency of aircraft wastewater testing to detect SARS-CoV-2 cases among travellers and to identify the viral genotype. It also highlights the low efficacy of the current control strategy for flights entering France from outside Europe, which combines a requirement to produce a vaccine pass and proof of a negative test before boarding.

Screening and Whole Genome Sequencing of SARS-CoV-2 Circulating During the First Three Waves of the COVID-19 Pandemic in Libreville and the Haut-Ogooué Province in Gabon.

Since the onset of the COVID-19 pandemic, the SARS-CoV-2 viral dynamics in Africa have been less documented than on other continents. In Gabon, a Central African country, a total number of 37,511 cases of COVID-19 and 281 deaths have been reported as of December 8, 2021. After the first COVID-19 case was reported on March 12, 2020, in the capital Libreville, the country experienced two successive waves. The first one, occurred in March 2020 to August 2020, and the second one in January 2021 to May 2021. The third wave began in September 2021 and ended in November 2021. In order to reduce the data gap regarding the dynamics of SARS-CoV-2 in Central Africa, we performed a retrospective genotyping study using 1,006 samples collected from COVID-19 patients in Gabon from 2020 to 2021. Using SARS-CoV-2 variant screening by Real-Time Quantitative Reverse Transcription PCR (qRT-PCR) and whole genome sequencing (WGS), we genotyped 809 SARS-CoV-2 samples through qRT-PCR and identified to generated 291 new genomes. It allowed us to describe specific mutations and changes in the SARS-CoV-2 variants in Gabon. The qRT-PCR screening of 809 positive samples from March 2020 to September 2021 showed that 119 SARS-CoV-2 samples (14.7%) were classified as VOC Alpha (Pangolin lineage B.1.1.7), one (0.1%) was a VOC Beta (B.1.351), and 198 (24.5 %) were VOC Delta (B.1.617.2), while 491 samples (60.7%) remained negative for the variants sought. The B1.1 variant was predominant during the first wave while the VOC Alpha dominated the second wave. The B1.617.2 Delta variant is currently the dominant variant of the third wave. Similarly, the analysis of the 291 genome sequences indicated that the dominant variant during the first wave was lineage B.1.1, while the dominant variants of the second wave were lineages B.1.1.7 (50.6%) and B.1.1.318 (36.4%). The third wave started with the circulation of the Delta variant (B.1.617). Finally, we compared these results to the SARS-CoV-2 sequences reported in other African, European, American and Asian countries. Sequences of Gabonese SARS-CoV-2 strains presented the highest similarities with those of France, Belgium and neighboring countries of Central Africa, as well as West Africa.